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Abstract 

Understanding and controlling the flow of heat is a major challenge in nanoelectronics. When a junction 
is driven out of equilibrium by light or the flow of electric charge, the vibrational and electronic degrees 
of freedom are, in general, no longer described by a single temperature ||l]46l. Moreover, characterizing the 
steady-state vibrational and electronic distributions in situ is extremely challenging. Here we show that 
surface-enhanced Raman emission may be used to determine the effective temperatures for both the vibra- 
tional modes and the flowing electrons in a biased metalhc nanoscale junction decorated with molecules Q. 
Molecular vibrations show mode-specific pumping by both optical excitationllSl and dc current||9l, with 
effective temperatures exceeding several hundred Kelvin. AntiStokes electronic Raman emission llTOl [TTI 
indicates electronic effective temperature also increase to as much as three times its no-current values at bias 
voltages of a few hundred mV. While the precise effective temperatures are model-dependent, the trends as 
a function of bias conditions are robust, and allow direct comparisons with theories of nanoscale heating. 
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When a bias, V, is applied to a nanoscale junction, the electronic distribution is driven out of 
equilibrium, with "hot" electrons injected from the source electrode into empty states eV above 
the Fermi level of the drain electrode. The electronic distribution rethermalizes on a scale set by 
inelastic electron-electron scattering. The Landauer-Biittiker approach assumes that dissipative 
processes are spatially removed from the junction region, so that the source and drain electronic 
distributions have identical temperatures [i 1 211 . On longer time and distance scales, electron-phonon 
scattering transfers energy to the lattice. If some current flows through molecules, inelastic pro- 
cesses can also transfer energy from the electrons to local molecular vibrational modes [|T3l [T4l . 
In the presence of optical radiation, the situation is even more complicated f lS] . Direct optical 
absorption and the decay of plasmon excitations produce hot electrons and holes in the metal with 
energies as large as those of the incident photons. Raman scattering processes, hot electrons, and a 
hot substrate can also excite molecular vibrational modes, even in the absence of direct bias-driven 
processes dH. 

Directly accessing the electronic and vibrational distributions in a nanoscale junction is ex- 
tremely difficult. Like others before us,[|3]-|5l|9l we assume that one may parametrize these (gen- 
erally complicated) distributions by effective temperatures. Previous experiments have used the 
stability of mechanical break junctions as a proxy for an effective ionic temperature ll5l [T6l [TTll . 

Raman spectroscopy is a more direct probe of energy distributions, and the effective tempera- 
ture of a particular vibrational mode, T^^, can be defined as[l9l [T8l[T9l 
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where and I^^ are the Stokes and antiStokes Raman intensities for that mode, is the incident 
laser frequency, hco^ is the Raman shift, and is a correction factor that accounts for the ratio 
of the antiStokes and Stokes cross-sections. In surface-enhanced Raman spectroscopy (SERS), 
plasmons enhance Raman response by orders of magnitude[10J. In this case, also accounts 
for differences in the Raman enhancement between the antiStokes and Stokes frequencies. We 
assume that A^, does not vary significantly with bias, since plasmonic properties are determined 
by metal geometry and dielectric function. Bias-driven changes in the intrinsic antiStokes/Stokes 
crossections would require Stark-like physics, and would depend strongly on molecular orientation 
and mode. 

Optical pumping may also elevate T^^ above the ambient temperature, T. If the Stokes process 
occurs faster than vibrational relaxation to the bulk substrate, the excited vibrational population 



will exceed the equilibrium average. Moreover, optically excited energetic electrons in the metal 
can couple to molecular vibrations. In this regime, I^^ would scale quadratically with incident 
laser intensity, while scales linearly as usualflH. 

Nanojunction devices open the possibility of observing bias-driven vibrational pumping via 
antiStokes/Stokes SERS measurements, as suggested theoretically [fT5ll and experimentally [|9l [T8l 
[T9l . Cross-sections for electron- vibrational pumping are expected to differ from mode to mode, 
as are the relaxation rates to the bulk substrate. Conversely, intra-molecular vibrational relaxation 
would push all vibrational modes toward a single effective temperature. 

Finally, in SERS systems with strong enhancements, a continuum of Raman scattering is ob- 
served from the conduction electrons in the metal [[TOl [TTl |20] - I22ll . Stokes scattering results from 
creation of electron-hole excitations, while antiStokes scattering requires excited electrons above 
the Fermi level of the metal. Crystal momentum kinematically limits the accessible Raman shifts, 
with disorder and the surface playing critical roles [[TTl |20l . There is no complete theoretical de- 
scription of the full electronic continuum[22]. However, the antiStokes Raman response at some 
shift e is expected [[Tn to be proportional to the joint density of states (JDOS) for electrons at 
energy e + E and holes at energy E. Taking the kinematic restriction and the band densities of 
states to be slowly varying over the antiStokes shift range of interest, and neglecting any resonance 
effects, we expect 



where the integral is over the conduction band, f{E, T) is the Fermi-Dirac function, and Tf^ is an 
effective temperature for the conduction electrons. Only Raman emission from the electrons that 
are proximate to the SERS "hotspot" is readily detected, since only that emission is appropriately 
enhanced. 

We perform Raman spectroscopy on pairs of gold nanojunction electrodes, with the fabrica- 
tion and measurement processes described in Methods. Prior to electromigration [[23|[ a monolayer 
of the molecule of interest is self-assembled on the electrodes. Two molecules were examined: 
a three-ring oligophenylene vinylene terminated in amine functional groups (0PV3, described 
further in Supplementary Information) and 1-dodecanethiol. Electromigration and measurements 
take place in vacuum with the substrate at 80 K. Such nanogaps are excellent SERS substrates [|24| . 
and have a long track record in tunnelling transport through single molecules [[141 [25l . In nano- 
junctions with measurable conductances, prior work[|7l|26l and the present experiments (see Sup- 
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plementary Information) demonstrate strong correlations in the time variation of the nanojunction 
conductance and "blinking" and spectral diffusion of the SERS spectrum. The tunnelling con- 
ductance probes a very small volume. Hence, correlations between the tunnelling conductance 
and the SERS signal imply single- or few-molecule SERS sensitivity. Final nanojunction conduc- 
tances range from 0.05 Go to 0.01 Go (where Go = 2e^ /h is the quantum of conductance) for 
junctions with molecules used in heating measurements, higher than that expected for an idealized 
gap bridged by such a molecule attached to both electrodes [|27l . For measurements of electronic 
temperature, nanojunction conductances ranged from 0.5 Go to 0.01 Go, with the highest conduc- 
tance nanoj unctions showing no signs of vibrational Raman indicating the absence of molecules 
in the nanojunction. For lower values of the conductance, direct metal-metal tunnelling continues 
to be responsible for much of the total conductance, but the temporal correlations between con- 
ductance and Raman spectra imply that at least some of the tunnelling current interacts with the 
molecule (or molecules) even if it (they) does not neatly bridge the gap between the electrodes 
(see Supplementary Information). An image of a typical nanojunction is shown in Fig. la, and the 
measurement scheme is shown in Fig. lb. 

Simultaneous optical and electrical measurements provide a wealth of information, as seen in 
Figure Ic, d. For this device we can observe optical and electrical pumping of individual vibrations 
of 0PV3. We also observe apparent heating of the electrons as evidenced by the rise in low 
wavenumber antiStokes continuum with increasing magnitude of bias voltage. For simplicity in 
presentation, we discuss the effective temperature analysis in terms of optical vibrational pumping, 
electrical vibrational pumping, and electron heating in turn. 

Optical pumping of vibrations {Tf^ > T at Voc = 0) was observed in nanojunctions for both 
test molecules. Figure 2a shows the evolution of the Raman response of a dodecanethiol junc- 
tion as a function of time. Several different discrete spectral configurations are observed, as seen 
on the Stokes side of the spectrum. For each of these configurations, we observe the same two 
modes (1099 cm^^ and 1475 cm^^) in the antiStokes spectrum, but with varying intensities. Ef- 
fective temperatures are inferred using Eq. (1), assuming the Raman enhancement is frequency 
independent and the cross-sections for Stokes and antiStokes scattering are equal {A^j = 1). Effec- 
tive temperatures range from 300 K to greater than 2000 K, demonstrating that vibrations may be 
driven far from equilibrium via optical excitation alone. Note that when the spectrum returns to a 
previously visited configuration, the effective temperatures also return to the previous values asso- 
ciated with that spectrum. The simplest explanation for the simultaneous changes in spectrum and 
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effective temperatures is that the molecule(s) in question undergo conformational changes relative 
to the Au electrodes. Such changes in conformation can affect SERS enhancement, vibrational 
relaxation times, electronic coupling between the molecule and electrodes, and possibly spectral 
diffusion as well. Any of these changes could result in different pumping rates. Note that the two 
modes presented here are not in equilibrium with one another, with an effective temperature differ- 
ence of more than 300 K between them. As shown in Fig. 2c, we observe a quadratic dependence 
of the antiStokes signal on incident power, as expected in the optical pumping regime. [[8]| 

Under DC bias we observe strong electrical pumping of specific vibrational modes in both 
0PV3 and dodecanethiol. Additional examples are given in Supplementary Information. Figures 
3a and d show the bias dependence of for different modes for two different devices with 
0PV3 assembled on them. We observe an approximately linear increase in T^^ with increasing 
\V\. This observation across multiple modes, multiple devices, and multiple molecules makes 
it very unlikely that this trend results from some bias-driven change in A,y. For the sample in 
Figure 3a-c we are signal limited below V = 0.2 V and observe no measureable antiStokes signal 
above our detector noise floor. There is a voltage asymmetry in conductance and T^^ between the 
positive and negative sweeps. On the positive V side, the 1625 cm^^ mode is only slightly higher 
in effective temperature than the 1317 cm^^ mode, while on the negative V side, the difference is 
considerably larger. This dichotomy is reproduced on multiple voltage sweeps on this device. One 
possible explanation for this sort of mode-to-mode variation and voltage asymmetry is that the 
molecule's electronic coupling is different for the two electrodes, so the effective cross-sections 
for electron-vibrational pumping depend on the direction of current flow. Such an explanation 
has been suggested for asymmetries observed in vibration-mediated inelastic electron tunneling 
spectroscopy! 28 1 . Figure 3b also shows another feature that we observe with some regularity 
in such junctions: a systematic shift in energy of about 15 cm^^ is observed for most (Stokes) 
vibrational modes as a function of V. This is an example of a Raman Stark shift [|29l . The decrease 
in overall Raman intensity observed on this particular positive bias sweep is attributed this to slow 
drift of the laser spot, and does not appear to affect the effective vibrational temperature since we 
are in a regime with no detectable optical pumping. 

We think it is unlikely that the vibrational heating results from heat flow from hot electrodes 
rather than the flow of charge through the molecule(s), for the following reason. If the electrodes 
acted as an effective high temperature bath, all molecular modes with any coupling to that bath 
should reach a similar steady state temperature. Instead, we observe only a few modes, and those 
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have disparate temperatures and heating rates. While plasmon dispersion in A,y is a concern, 
numerical modeling of such j unctions [|24ll shows trends (stronger resonances to the red; resonances 
hundreds of meV broad in energy) inconsistent with that as a dominant effect. While we cannot 
rule out a contribution due to energy flow from the electrodes to the molecule(s), we believe 
pumping of vibrations due to current[28J is more consistent with these observations. 

For the device shown in Figure 3d, e, and f, both electrical and optical vibrational pumping are 
observed. The clear 1815 cm^ antiStokes mode has measurable antiStokes signal at zero bias, 
with a temperature more than 300 K above the substrate temperature. Until V exceeds hcu^, no 
increase in T^^ is observed, as expected. Once V exceeds this threshold, once again an approxi- 
mately linear increase in T^^ is observed with increasing V. These observations are consistent with 
theoretical treatments of electron- vibrational scattering during off-resonant conduction[|4l|30l. We 
do not observe any signs of vibrational cooling at high biasesjH [H; nor do we see T^^ ~ \/V 
at lower biases, as would be expected for cooling of the local modes via relaxation to bulk 3d 
phonons[3|. This suggests that local vibrational relaxation to the substrate takes place via other 
means, such as coupling to phonons of reduced dimensionahty. We again observe some asym- 
metry in conductance (Figure Id) and T^;^ between the positive and negative sweeps. The vibra- 
tions at 1480 cm^^ and 562 cm^^ show a strong asymmetry similar to the conductance, while the 
1815 cm^^ mode shows very little asymmetry. This again suggests that effective cross-sections 
for electron- vibrational pumping differ from mode to mode. 

In Figure 4 we examine heating of the electrons as a function of V. Effective temperatures are 
found by fitting the antiStokes continuum spectra to Equation 2. AtV = Vmax, two fit parameters 
are used, T^^ and an overall amplitude. At other bias voltages, the amplitude is then fixed, leaving 
Tg^ as the only free parameter to fit at each V. Representative spectra and fits are shown in Figs. 
4c, f. Error estimates for T^^ correspond to a doubling of from the best fit. In Figure 4a we see 
a roughly 150 K change in temperature over a 400 mV change in bias. There is asymmetry in both 
the temperature and current as a function of V. One possible explanation for this would be that 
only one electrode gets hot at a time, depending on the direction of current flow, and the thermal 
coupling to the bulk is different for each side of the nanogap. An increase of T^^ roughly linear in 
V is consistent with a general argument based on steady state heat flow in the electron system||3l. 

The observation of effective electronic heating with bias in nanoscale junctions is important. 
While it is expected on general grounds flSj, the Landauer-Biittiker approach to conduction, com- 
monly used to model such junctions, does not account for such an effect. Additional experimental 
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consequences would be expected, including broadening of inelastic electron tunneling spectra at 
high biases under similar current densities, and signatures of heating in the junction current noise. 
Investigations in these directions are ongoing. 

There are systematic uncertainties in the effective vibrational and electronic temperatures (due 
to the unknown values of and the lack of a detailed theory of the electronic Raman contin- 
uum). However, the trends captured in the data (the dependence of the effective temperatures on 
the bias conditions, for example) are robust. These experiments demonstrate that it is possible 
to extract highly local information about the effective ionic (vibrational) and electronic tempera- 
tures within nanoscale junctions under electrical bias. The resulting data are among the first direct 
measurements that will allow specific tests of theoretical treatments [[B [21 |4l|6l of thermal dissipa- 
tion at the nanoscale. It is clear that there can be great variation in electron-vibrational couplings 
and vibration-bulk relaxation rates; with weak relaxation to the bulk, electrical pumping could 
access vibrational temperatures comparable to eV/kB- These experiments may open the door to 
deliberate control and engineering of such effects. 
Methods 

Samples were fabricated on N-type Si with 200 nm of thermal oxide. Nanowires approximately 
120 nm wide were defined by e-beam lithography and e-beam evaporation of 0.5 nm Ti and 15 nm 
of Au. The nanowires were approximately 600 nm long and connected to larger electrical leads and 
contact pads. After fabrication samples were O2 plasma cleaned and soaked in either 0.1 mg/mL 
0PV3/THF or 10 mM 1-dodecanethiol in ethanol for 24 hours to form a self-assembled monolayer 
on the Au surface. After assembly samples were wire bonded to a ceramic chip carrier and placed 
in an optical cryostat. Samples were then cooled to 80 K in vacuum. 

Nanogaps were formed in the nanowires via the electromigration process which has been 
studied extensively. Gap formation was performed with an automated computer program which 
stopped the electromigration process once the electrical conduction dropped below 2e^//i = Go - 
the quanta of conductance. 

After electromigration optical measurements were performed by illuminating the sample sur- 
face at 785 nm with a Gaussian laser beam with a spot size of 1.9 /xm, incident intensity of 
23 kW/cm^, and polarization aligned along the nanowire. The laser beam was rastered over the 
sample and a spatial map of the Raman response of the Si substrate was measured at 520 cm~^ 
as described in Ref. i7j. The Au electrodes of the nanogap attenuate the silicon signal allowing 
precise positioning of the laser over the nanogap. The SERS signal of the nanogap was then mea- 
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sured. The SERS signal may also be spatial mapped to confirm that it originates from the nanogap. 
The spectra were all corrected to adjust for the spectral sensitivity of the system over the measured 
wavelength range, by calibrating to a known blackbody source (T=2960 K). Integration times for 
SERS spectra during voltage sweeps ranged from 5 s to 10 s depending on signal strength. 

Electrical measurements of the nanogap were performed using a current amplifier and lock-in 
amplifier as illustrated in Figure lb. The lock- in amplifier applied a 10 mV AC signal to the sample 
and measures the first harmonic response (oc dl/dV). A summing amplifier is used to add a DC 
signal (V) to the lock-in signal allowing the current amplifier to measure /. In this configuration 
both / and dl /dV can be measured simultaneously as functions of V . 

To avoid concerns about thermal broadening or extrinsic (detector-limited) broadening of Ra- 
man peaks, the effective temperature analyses were performed based on absolute peak height 
above baseline rather than integrated intensity over some wavenumber range. Error estimates take 
into consideration the fact that the peak heights have a noise, where is the number of CCD 
counts over baseline. 

[Supplementary Information] accompanies this paper at www.nature.com/naturenanotechnology. 
Reprints and permission information is available online at http://npg.nature.com/reprintsandpermis sions/ . 
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FIG. 1. Measurement overview, a) SEM false color image of a typical nanogap. b) Schematic of electrical 
and optical measurement, c) Waterfall plot of Raman response (in CCD counts) of an 0PV3 junction as a 
function of DC bias, V. The antiStokes spectrum shows strong dependence on V while the Stokes side is 
relatively constant as a function of bias. The strong Stokes peak at 520 cm^^ is from the silicon substrate, 
d) Simultaneously measured current and differential conductance as a function of V. 
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FIG. 2. Optically driven vibrational pumping, a) Raman response of a dodecanethiol-coated junction 
as a function of time under zero bias. Scale bars are in thousands of CCD counts. The junction switches 
stochastically between several stable configurations, each with characteristic spectra that exhibit strong 
optical pumping of vibrational modes at 1099 cm~^ (C-C stretch) and 1475 cm^^ (CH2 or CH3 deform), 
b) Spectra from given time slices in (a). We note that the higher energy peak spectrally diffuses between 
1475 cm^^ and 1525 cm~^ and the lower energy peak between 1099 cm^^ and 1110 cm^^. The temperature 
of the 1099 cm~^ mode in Kelvin is 769, 1722, 532, 754, 512, 775 with increasing time. The temperature 
of the 1475 cm"^ mode in is 701, 815, 658, 729, 648, 746. We note that at 70 s a vibration at 1450 cm"^ is 
observed at 2161 K. All temperatures have an uncertainty of it 10 K. c) Power dependence of Stokes (black) 
and antiStokes (blue) signal taken from a different device also showing optical pumping. As expected, in 
the optical pumping regime the Stokes signal is linear (solid line) in laser power while the antiStokes signal 
increases quadratically (solid line) in laser power. Error bars indicate uncertainty in signal due to read and 
shot noise in CCD. 



10 




-0.5 -0.25 0.25 0.5 .igoo .1400 1400 1600 

V Raman Shift [cm'''] 



700 



600 



0) 500 



400 




200 
180 


1 


0.4- 
0.3 - 


160 
140 




0.2 - 


120 




0.1 - 


100 




vo- 


80 




-0.1 - 


60 
40 




-0.2- 


20 


1 


-0.3- 







-0.4- 



10000 

9000 

8000 

7000 " 
'uj 

6000 I 

c 
ro 

5000 E 
ro 

d: 

4000 
3000 



500 mV rt 




400 mV A 


^ 


300 mV 




200 mV 




lOOmV 




OmV 






-1600 -1400 1400 1600 

Raman Shift [crn^] 



ij^ — ^ 


400 mV . 




320 mV 1 




240 mV I 




160 mV 


A — — 


80 mV 


»--.. — - 


OmV 



-1800-1300 -800 800 1300 1800 
Raman Shift [cm"^] 



-1800-1300-800 800 1300 1800 
Raman Shift Icm^l 



FIG. 3. Electrically driven vibrational pumping, a) Effective vibrational temperature as a function of V 
for two 0PV3 modes: 1317 cm^^(red) and 1625 cm~^(blue). Error bars indicate uncertainty in temperature 
due to antiStokes amplitude measurements. Inset) IV curve for this device, b) Raman response of this 
device as a function of V. c) Sample spectra for given voltage. All Stokes and antiStokes are plotted on 
the same scale. Full amplitude corresponds to 235 (10,000) counts for antiStokes (Stokes), d) Effective 
vibrational temperature as a function of V for three 0PV3 modes: 1815 cm' 'Hblack), 1480 cm'^red) and 
562 cm^^(blue). Error bars indicate uncertainty in temperature due to antiStokes amplitude measurements. 
Simultaneous optical and electrical vibrational pumping are observed for the 1810 cm^^ mode as Teff is 
much greater than 80 K at F = 0. e) Raman response of this devices as a function of V. Blue indicates 
counts and red indicates 200 (5000) counts for antiStokes (Stokes), f) Sample spectra for given voltage. All 
Stokes and antiStokes are plotted on the same scale. Full amplitude corresponds to 350 (10,000) counts for 
antiStokes (Stokes). 1 1 




FIG. 4. Electronic heating under bias, a&d) Effective temperature (blue) and dissipated electrical power 
(red) for two different devices. The device in (a, b, c) shows no molecular Raman peaks, and is considered 
a "clean" junction, while (d, e, f) is the same 0PV3 device that was analyzed in Figs. 1 and 3. Insets: I — V 
curves for these devices. Error bars are described in text, b&e) Raman response for these devices, c&f) 
Sample spectra and best fits given by Equation 2 for a given voltage. 
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Supporting Information: 
Vibrational and electronic heating in nanoscale junctions 

I. CORRELATED ELECTRICAL CONDUCTANCE AND RAMAN RESPONSE DATA 

Here are several examples of correlated electrical conductance and Raman response from nano- 
junctions with OPV3 assembled on them. In the first three examples (Figures SI -3) the electrical 
conductance is on the order of 0.01 Go with conductance fluctuations on the order of 0.001 Go- 
The conductance fluctuations are correlated with large changes in the Raman response measured at 
the nanojunction. This is consistent with a single molecule switching between different conforma- 
tions in the nanojunction. In all of these cases, we expect that the molecule is not neatly bridging 
the gap, but does play a role in total electrical conduction. It has long been established[S2-S4] 
that conduction in such nanojunctions takes place via tunneling. Because of the exponential de- 
pendence of tunneling on distance, the dominant volume for current flow is of molecular scale. As 
argued previously [S5], coincident conductance and Raman changes imply (1) that at least some 
portion of a Raman-active molecule is influencing the tunneling electrons; and (2) that Raman- 
active molecule is a significant (in many cases dominant) contributor to the total Raman signal. 

At higher nanojunction conductances (0.5 Gq, such as those observed in Figure S4), we ob- 
serve similar correlations between conduction and Raman response, but the conductance fluctua- 
tion magnitude is much larger, as one might expect. At this conductance level, conduction is likely 
dominated by an atomic-scale metallic junction with transmission less than one. (One example of 
such a junction would be the beginnings of a tunneling contact, with two tip atoms slightly farther 
apart than their equiUbrium lattice spacing.) Again, the conduction path is highly localized, with 
a transverse dimension of molecular or atomic dimensions. Still, given the correlation between 
Raman response and conductance, at least one Raman-active molecule must be significantly cou- 
pled to this current path. We note that it is very unlikely that the unusually high conductance of 
this particular junction results from many molecules in parallel. If this was the case, the observed, 
correlated conduction and Raman response would imply many molecules changing their configu- 
rations simultaneously, which is very unlikely. One would expect conductance fluctuations due to 
individual molecular motions to be much smaller such as those seen in nanojunctions with con- 
ductances of 0.01 Go; however this is not the case. Instead, high conductance nanojunctions are 
likely bridged by or strongly electronically coupled to just a part of a molecule, as the gaps are too 
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small (based upon the conductance) to fit an entire molecule. 



n. JUNCTION CONFIGURATIONS AND CONTRIBUTIONS TO THE CONDUCTANCE 



Figure S5 presents schematic examples of possible nanojunction configurations. Fig. S5i 
shows the idealized single-molecule junction that has been considered for more than fifteen years, 
with a single molecule neatly bridging and bound to both ends of a nanoscale interelectrode gap. 
In such a geometry, it is clear that interelectrode conduction would have a dominant contribution 
from current flow through the molecule. 

However, as mentioned in the main text, this ideal configuration is unlikely in electromigrated 
junctions, since the electromigration process (unlike mechanical break junctions or scanning tun- 
neling microscope junctions, for example) does not control interelectrode distance with atomic 
precision. Rather, the closest interelectrode distance is determined by details of individual junc- 
tions and the breaking procedure. 



Fig. S5 3-d are examples of other nanojunction geometries, in which the closest interelectrode 



distance is not the length of a molecule of interest. In Fig. S5 f and d, only a portion of the 



molecule of interest may be located at the point of closest interelectrode separation, while steric 



effects prevent the molecule from forming strong bonds with both electrodes. In Fig. S5 the 
molecule is bound to both electrodes, but a closer point of interelectrode separation exists, and due 
to the exponential decay of tunneling with distance, the total conductance would have a dominant 
contribution from direct metal-metal tunneling. These configurations, while only schematic, show 
clearly that many reasonable junction arrangements are possible that can have total conductances 



higher than the idealized situation in Fig. S5i, but nonetheless involve contributions from current 
flow through a molecule. An example from the literature of such a situation is thought to be 
presented in Tal et a/.[S6]. In that work, a benzene molecule is thought to lie down flat between 
two Pt electrodes in a mechanical break junction, leading to high conductance (comparable to 1 
Go), but still with molecular influence on the current. 



III. ADDITIONAL OPTICAL VIBRATIONAL PUMPING DATA 

An example of optical vibrational pumping in 0PV3 is shown in Figure S6. We observe several 
discrete spectral configurations. Each configuration has varying levels of optical pumping, with 
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the most noticable change occuring at 150 s, when several lower energy vibrations appear in the 
antiStokes signal. 

IV. ADDITIONAL ELECTRICAL VIBRATIONAL PUMPING DATA 

An additional example of electrical vibration pumping in OPV3 is shown in Figure S7. The 
two modes observed 880 cm~^ and 410 cm~^ have temperature differences of almost 150 K at 
V — 0.4 V. This device also exhibited optical pumping of both modes resulting in zero bias 
temperatures of 220 K and 120 K respectively. As expected both modes to not show any electrical 
pumping until V exceeds the vibrational energy. At 280 mV a sharp change in conductance is 
observed and no pumping is observed. The conductance changes again at 320 mV and pumping is 
observed but with less AS intensity than previously. This is strong evidence the importance of the 
local environment to pumping cross-sections. 

V. ADDITIONAL ELECTRON HEATING DATA 

An additional example of electrical heating is shown in Figure S8. This sample had more sym- 
metric I—V curves and power dissipation than other samples presented. The effective temperature 
is observed to be more symmetric as well when considering the size of the error bars. 

VI. RAMAN STARK EFFECT 

Figure S9 provides a closer look at the Raman Stark effect seen in Figure SB of the main text. 
A clear shift of 13 cm~^ can be observed between zero bias and 0.5 V. All vibrational modes 
present experience similar shifts. These shifts are reproducible on subsequent voltage sweeps of 
this junction. 

VII. OPV3 SYNTHESIS 

Tetrahydrofuran (THF) was distilled from sodium benzophenone ketyl. Triethylamine (TEA) 
was distilled from CaH2 under N2. Triethylphosphite (98%, Aldrich), 4-aminostyrene (90%, 
Aldrich), and palladium(II) acetate (98% Aldrich) were used as received. SiUca gel plates were 
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250 mm thick, 40 F254 grade obtained from EM Science. Silica gel was grade 60 (230 - 400 
mesh) from EM Science. 

Anhydrous DMF (120 mL) was placed in a 250 mL round bottom flask and two freeze- 
pump-thaw cycles were performed to ensure the removal of oxygen. A large screw-cap tube 
was charged with 2,5-didecoxy-l,4-diiodobenzene (10.4 g, 16.2 mmol),[Sl] tetrabutylammonium 
bromide (15.7 g, 48.6 mmol), Pd(0Ac)2 (0.40 g, 1.62 mmol), and K2CO3 (4.48 g, 32.4 mmol). 
4-Aminostyrene (5.00 g, 32.4 mmol) was added to the DMF and the resulting solution was can- 
nulated into the screw-cap tube. The tube was sealed with the Teflon cap and the tube was heated 
to 100 °C for 1 d. After cooling, the crude reaction mixture was poured into water, extracted with 
Et20, dried over MgS04, and then the solvent was removed under reduced pressure. The crude 
product was then purified by silica gel chromatography using 1:1 CH2CI2 : hexanes as eluent. The 
solid was then recrystallized from CH2CI2 : hexanes to yield 1.57 g (16%) of an orange solid. IR: 
3729, 3470, 3375, 3 207, 3025, 2925, 2851, 2648, 2321, 1976, 1622, 1599, 1516, 1466, 1413, 
1254, 1172, 962, 895, 812 cm^^ NMR (400 MHz, CD2CI2) 5 7.26 (d, J = 8.3 Hz, 4H), 7.18 
(d, J = 16.5 Hz, 2 H), 7.01 (2 H), 6.96 (d, J = 16.5 Hz, 2 H), 3.95 (t, J = 6.5 Hz, 4 H), 3.73 (br 
s, 4 H), 1.78 (m, 4 H), 1.45 (m, 4 H), 1.34 (br m, 24 H), 0.80 (m, 6 H); (100 MHz, CD2CI2) 
5 151.35, 147.02, 128.95, 128.91, 128.18, 127.21, 120.02, 115.51, 110.58, 70.10, 32.51, 30.26, 
30.19, 30.14, 30.06, 29.96, 26.88, 23.29, 14.48. HRMS calcd for C42H60N2O2: 624.4655, found: 
624.4662. 



The structure of the 0PV3 Molecule is shown in Figure SIO 
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FIG. S 1 . Waterfall plot of Raman spectrum (Is integrations) and electrical conduction measurement for 
an 0PV3 devices. The Raman response is observed to change whenever a change in conductance occurs. 
Colorbar indicates Raman intensity in CCD counts. 
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FIG. S2. Waterfall plot of Raman spectrum (Is integrations) and electrical conduction measurement for 
an 0PV3 devices. The Raman response is observed to change whenever a change in conductance occurs. 
The number of Raman vibrational modes observed also changes with conductance level. Colorbar indicates 
Raman intensity in CCD counts. 
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FIG. S3. Waterfall plot of Raman spectrum (5 s integrations) and electrical conduction measurement for 
an 0PV3 devices. The Raman response is observed to change whenever a change in conductance occurs. 
Colorbar indicates Raman intensity in CCD counts. 
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FIG. S4. Waterfall plot of Raman spectrum (Is integrations) and electrical conduction measurement for 
an 0PV3 devices. The Raman response is observed to change whenever a change in conductance occurs. 
Colorbar indicates Raman intensity in CCD counts. 
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FIG. S5. Cartoon of possible junction configurations, using the 0PV3 molecule as an example (saturated 
side chains omitted for clarity), (a) The idealized single-molecule junction, a configuration unlikely in 
these experiments since interelectrode separation is not precisely controlled, (b-d) Alternative junction 
configurations, in which interelectrode conduction would be expected to include a contribution from current 
interacting with the molecule, as well as direct metal-metal tunneling. 
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FIG. S6. Raman response of an 0PV3 junction as a function of time under zero bias. Blue indicates 
counts and red indicates 100 (8000) counts for antiStokes (Stokes) sides. Integration time is 1 s. The 
junction switches stochastically between several stable configurations, each with characteristic spectra that 
exhibit strong optical pumping of different vibrational modes. 
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FIG. S7. A) Effective vibrational temperature as a function of V for two 0PV3 modes: 880 cm~^(red) and 
410 cm~^(blue). Error bars indicate the uncertainty in inferred effective temperature due to the statistical 
limitations of the antiStokes amplitude measurements. Inset) IV curve for this device. B) Raman response 
of this device as a function of V. Blue indicates 10 (2500) counts and red indicates 250 (10,000) counts for 
antiStokes (Stokes). The strong Stokes peak at 520 cm~^ is from the Si substrate. C) Sample spectra for 
given voltage. All antiStokes (Stokes) spectra are plotted on the same scale. Full amplitude corresponds to 
1200 (18,000) counts for antiStokes (Stokes). 
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FIG. S8. A)Effective electronic temperature(blue) and dissipated electrical power (red) for a device with 
very little vibrational Raman activity. Error bars are described in text. Inset) I — V curves for these devices. 
B) Raman response for these devices. Blue indicates counts and red indicates 350 counts. C) Sample 
spectra (blue) and best fit given by Equation 2 in main text (green) for a given voltage. 
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FIG. S9. A) Raman response of the same device shown in Figure 3a in main text. A shift of about 15 cm~ 
is present for many spectral lines. Blue indicates 2500 counts and red indicates 10000 counts. B) Zoom in 
of Raman response for the mode centered at 1523 cm^^.Blue indicates 2500 counts and red indicates 7000 
counts. C) Sample spectra for given voltage. All spectra are plotted on the same scale with a base line of 
3000 counts substracted. Full amplitude corresponds to 4000 counts. The peak at 1523 cm^^ can clearly 
be seen systematically shifting to lower energy at higher voltages reaching 1510 cm^^ at 500 mV. 
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FIG. SIO. 0PV3 Molecule 
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